Statistics dependence of the entanglement entropy 
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The entanglement entropy of a distinguished region of a quantum many-body system reflects the entanglement 
present in its pure ground state. In this work, we establish scaling laws for this entanglement for critical quasi- 
free fermionic and bosonic lattice systems, without resorting to numerical means. We consider the geometrical 
setting of .D-dimensional half-spaces which allows us to exploit a connection to the one-dimensional case. 
Intriguingly, we find a difference in the scaling properties depending on whether the system is bosonic — where 
an area-law is first proven to hold — or fermionic, extending previous findings for cubic regions. For bosonic 
systems with nearest neighbor interaction we prove the conjectured area-law by computing the logarithmic 
negativity analytically. We identify a length scale associated with entanglement, different from the correlation 
length. For fermions we determine the logarithmic correction to the area-law, which depends on the topology 
of the Fermi surface. We find that Lifshitz quantum phase transitions are accompanied with a non-analyticity in 
the prefactor of the leading order term. 



The occurrence of critical points at zero temperature holds 
the key to the understanding of several phenomena in quan- 
tum many -body systems in the condensed matter context (H. 
Quantum criticality is accompanied by a divergence of the 
typical length scale, the correlation length. These long-range 
correlations come along with genuine entanglement in the 
ground state, grasped by the entanglement entropy Eg = 
S(tr\j[[p\). This is the entropy of the reduced density ma- 
trix that is obtained when tracing out the degrees of freedom 
outside a distinguished region A, quantifying the degree of 
entanglement between this region and the rest [2-19]. 

This notion of the entanglement or geometric entropy — and 
more to the point its scaling behavior abstracting from details 
of the model — has enjoyed a strong revival of interest recently, 
partially driven by intuition from quantum information theory: 
previously conjectured scaling laws in higher dimensions H, 
relating the entanglement entropy to the boundary area — not 
the volume — of the region, have been rigorously established 
using quantum information ideas [4-6]. This was followed by 
observations of violations of such area-laws 0]. The entan- 
glement entropy has in its non-leading-order behavior inter- 
estingly been linked to the topology of the system Jit], using 
ideas of topological quantum field theory, and been studied 
under time evolution [9]. Partly, this renewed interest is trig- 
gered by the implications on the simulatability of quantum 
systems using density-matrix renormalization approaches: the 
entanglement entropy quantifies in a sense the relevant num- 
ber of degrees of freedom to be considered ll3i [Toll . 

If entanglement is to reflect critical or non-critical prop- 
erties of quantum many-body systems, an area-relationship 
might of course be expected to hold or not, depending on 
whether the two-point correlation functions diverge. One 
might be tempted to think that entanglement could yet be seen 
as an indicator of criticality in the same sense. Intriguingly, it 
turns out that the situation is more complex. As we will also 
see, even for critical systems, an area-relationship can hold, 
despite a divergent correlation length (as can also be observed 
in projected entangled pair states, satisfying an area-law by 
construction ifTTll '). In this work, we demonstrate with a fully 
analytical argument that it can depend on the statistics of the 




FIG. 1: Spectra of critical nearest-neighbour Hamiltonians (top: 
bosonic, bottom: fermionic) in one and two dimensions. The spec- 
trum of the individual decoupled chains is given by the spectrum of 
the full Hamiltonian along the first coordinate <pi. The topology of 
the set of solutions to \ v — — the Fermi-surface in the fermionic 
case — provides an intuition as to why the scaling behavior of entan- 
glement is different for fermions then for bosons in D > 1, see text. 



system — bosonic or fermionic — whether an area-relationship 
holds or is in fact violated. In this way, we resolve the key 
open question: "What happens in the critical bosonic case?" 
In the process we confirm some conjectures based on numeri- 
cal findings for small system size ll2T fl2i[T3ll and refute others, 
such as the conjecture of a break-down of an area law for crit- 
ical bosons in D > 1 in Ref. |fl4ll . 

Here, we establish first analytical scaling laws for critical 
bosonic systems. We achieve these results for the geometrical 
setting of a half-space in D-dimensions, completing the pro- 
gram initiated in Ref. 15]]. These findings are compared with 
fermionic half-spaces, complementing recent results on cubic 
regions in Refs. |7f|, a nd in a fashion consistent with numer- 
ical work in Ref. Il2l [l3ll . We treat bosonic and fermionic 
systems on the same footing - in terms of Majorana operators 
for fermions and canonical coordinates for bosons. We hence 
provide a unified and complete framework for entanglement 
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scaling in critical quasi-free systems with that geometry. 

The setting. - We consider cubic lattices of spatial dimen- 
sion D, C = [1, . . . , N] xD , \C\ = N D , equipped with peri- 
odic boundary conditions and study ground states q of Hamil- 
tonians that are quadratic forms of either bosonic or fermionic 
operators. The geometric setting is that of a half-space, distin- 
guishing w.l.o.g. the first spatial direction, considering a sub- 
system A = [1, . .., M] x £', C = [1, . . ., N}* ' 1 , and its 
entanglement with the rest B = C\A = [M + l, . . . , N] x 
We can hence make use of an idea of exploiting transverse 
momenta. This geometrical setting has notably numerically 
been assessed with respect to local spectra and simulatability 
issues in the seminal work Ref. yfl. 

When we say that (i) this entanglement as quantified by the 
von-Neumann entropy Es(M,N) = S(qa) of the reduced 
state Qj\ = trg[g] satisfies an area-law, we mean that for M = 
N/2 the entanglement entropy satisfies 

£ := lim E S {M,N)/N D - 1 < const, 

N— >oo 

i.e., it scales at most like the boundary area of A. 

For systems violating the area-law, £ = oo, we will (ii) 
study the exact form of the encountered logarithmic diver- 
gence in M, 

£ = const x log M + o(logM). 

We will subsequently discuss physical systems that are de- 
scribed by Hamiltonians of the type 



H = 



where operators di are either bosonic or fermionic and vec- 
tors i = - ■ - id) € C label individual sites of the cubic 
lattice. To ensure hermiticity we demand the real coefficients 



to satisfy Ai j = Bi j = 
bosons, and Aij = —Bi 



Aj i and Ci. 
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for fermions. Furthermore, we assume translational invari- 
ance and periodic boundary conditions (all coupling matrices 
depend only on the difference i — j and are cyclic matrices). 
We will lead the discussion in terms of hermitian operators 
r = (ii, i|£|,pi, ... ,p|£|)* (in a mild abuse of nota- 
tion, the transposition refers to the tuple, not to operators) de- 
fined by £i = (di + d\)/y/2 and pi = —i(di — d\)/\/2. In 
the bosonic case they are indeed position and momentum op- 
erators fulfilling the canonical commutation relations (CCR). 
In turn, for fermionic operators, they are so called Majorana 
operators fulfilling the canonical anti-commutation relations 
(CAR). We will assume isotropic couplings for fermions, 
C = 0, and coupling only in position for bosons, A — C — 1. 
In order not to obscure our main point, we will not con- 
sider the straightforward but cumbersome generalization to 
anisotropic or momentum couplings. The Hamiltonian now 
reads (V B := A + C, V F := A) 



Hr 
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V B 
1 



V F 

-vi o 



r (1) 



for bosons and fermions, respectively. Whenever we may treat 
both species equally, we denote by V the coupling in position 
V = Vb for bosons and V = Vf for fermions. 

Starting from the Hamiltonians above, their spectrum and 
respective ground states are found in the usual way by diago- 
nalizing Hb through symplectic transformations, transforma- 
tions respecting the CCR, and Hp by orthogonal transforma- 
tions, transformations respecting the CAR. As matrices V are 
cyclic, the bosonic spectrum is given by and the fermionic 
spectrum by |Afc|, where = X^e£ ^ cos(27rfei/n), k £ C. 

The half-space geometry allows for a transformation of 
both Hamiltonians to a system of mutually uncoupled one- 
dimensional chains while respecting the CCR and CAR, but 
notably, changing the local properties of the systems forming 
the individual chains. To this end consider the transformation 
r = (0@0)q to a new set of operators q, where the \C\ x \C\ 
matrix O is given by Oij — S^^Oi' j' . Here and in the fol- 
lowing we write vectors % = (if ■ ■ id) £ £ as i = (i\i'), 
i! = (12 ■ ■ - id) £ Now, the transformation O acts on the 
first coordinate as the identity and is thus local with respect 
to the bipartition A\B, i.e., it does not change entanglement 
properties. In order to respect the CAR and CCR the matrix 
O needs to be orthogonal. Then, the Hamiltonians H B / f rea d 
in coordinates q just as in Eq. (Q~|i with modified coupling ma- 
trices V i— > O l VO. As cyclic matrices commute and may be 
diagonalized by the same Fourier transformation, we have the 
explicit form 



k'ec 



where the A 



N x 27rifeiZi/iV 

771 , 

1 fci=i 



N 



tkik') = ^k, k G C, are the eigenvalues of V. 
Now, define the \C'\ x \C'\ matrices V(h) as those matri- 
ces obtained from V by keeping the first coordinate fixed: 
(V(li))i> = Vn^iy Then the Afc'(Zi) are the eigenvalues of 
the V(li), which are all cyclic and can thus all be diagonal- 
ized by the same orthogonal matrix. Choosing O to be this 
matrix yields (O t VO)^ lk '- ) ^ :jlp , ) = 6 k ', p ,\ k <(ii -ji), ^.mo- 
mentum space representation of the coupling in all but the first 
coordinate. In this representation the Hamiltonian is a sum of 
mutually uncoupled one-dimensional chains labeled by 
k'. Each chain is described by a Hamiltonian of the form 
as in Eq. (fl~|i with N x N cyclic coupling matrices V(k ), 
(V(k')) h = A fc <ii). We will write (V(<p'))^ = \ v (ipx), 
4>i = 2-kIi/N, tp' d = 2nk' d /N, for the infinite system. After 
this decoupling procedure the entanglement between A and B 
is now given by a sum of the entanglement between the sites 
[1, . . . , M] and [M + 1, . . . , N] of the individual chains: 



8= limV^ 

k' 



E s (ip') 



dtp' 



[0,2tt] 



(27T) 
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This will be the starting point for the following discussion. 

Fermions. - We start with investigating case (ii) above: The 
asymptotic behavior in M after taking the limit TV — > oo l20ll . 
frequently referred to as the double scaling limit lfl7ll . For 
each chain ip' we now need to compute the entanglement be- 
tween the first M sites and the rest of the chain (in real space, 
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not momentum space). The asymptotic behavior in M of 
this entanglement can be obtained from the so-called symbol 
g v {<Px) = sgn(A v <<?i)) of the chain Old]: Each E s (tp') 
is determined from the continuity properties of g v i as a func- 
tion of ipi. For fixed ip' it corresponds to a one-dimensional 
isotropic fermionic model, for which the asymptotic form of 
the entanglement has been obtained in Refs. fun . 



Es(<p') = 



sty) 
6 



\oz 2 M + c{ip') + o{\ogM), (2) 



where the function c(ip') does not depend on M and the inte- 
ger s(cp') is the number of discontinuities of g v > as function 
of tpi in the interval [0, 2ir). All chains with s(cp') > are 
critical as finding discontinuities in the symbol is equivalent 
to having a vanishing energy gap above the ground state. This 
means that, depending on the Fermi-surface (the set of solu- 
tions to \ v = 0, see Fig.Q]), one finds a continuum of chains 
ip' that are critical. This situation is in contrast to the situation 
encountered when considering bosonic systems as we will see 
below. From (f2]), we find the asymptotic behavior in M as 



log 2 M^ av ^ a ) 
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^g) +0 (logM), (3) 



where we defined v($a) = f$d.<p' as the volume of the set 
= {if 1 : s(ip') = a}, so the set associated with exactly 
a discontinuities. Hence, we do encounter a logarithmic di- 
vergence in M of the entanglement entropy and the pre-factor 
depends on the topology of the Fermi-surface: g v i exhibits 
discontinuities at points where \ v = 0, i.e., on the Fermi- 
surface. If the Fermi surface is of measure zero (i.e., the set 
of solutions to X v = is countable, as for example in the 
critical bosonic case, see Fig. [TJ, we have v($> a ) = and the 
system obeys the area law, £ = const. Consider as an ex- 
ample the case of a nearest neighbor Hamiltonian with cou- 



pling Vij 



aS, 



diSt(i,i),l: 



in which case the symbol 



corresponds for fixed ip' to that of the isotropic XY model 
with transverse magnetic field h((p') = 1 + 2a J2 d C0S ( L P'd)- 
For this model, the non-leading order term was obtained em- 
ploying Fisher-Hartwig type methods Oil- It reads c(<p') = 
log 2 (l — h 2 ((p')/A)/6 + cq and vanishes if the chain <p' is 
non-critical. The constant cq is independent of the system pa- 
rameters. The number of discontinuities is s((p') = 2 for 

<P' G *2 = {<p' £ [O^)^ 1 : |-L + ^cos(^)| < l} 



and zero otherwise, i.e., the sum over cr in ((3J consists only of 
the <T = 2 term as all others are zero. Thus 



£ 



3(2tt) d ~ 1 



log M 



(27T) 
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■o(logM), 



where for D = 2 and the critical case \a\ > 1/4, we find 
^($2) = 2arccos (l/(2|a|) — l), i.e., the prefactor depends 
on the coupling parameter a. For non-critical models in the 
isotropic setting at hand the set $2 is empty and there is no 
entanglement. There is no universal non-leading order term 



as proposed in Ref. [8[] related to the conformal charge, due to 
the specific geometric setting of a half-space. 

At this point, it is interesting to discuss the behavior of the 
entanglement entropy under Lifshitz phase transitions. They 
are topological quantum phase transitions of fermionic sys- 
tems due to a change of the topology of the Fermi surface, oc- 
curring for example in d-wave superconductors 0L see a lso 
Ref. 11311 . The previous considerations immediately allow us 
to argue that a Lifshitz transition accompanied by a change of 
topology of the Fermi surface is reflected by a non-analyticity 
in the prefactor of the entanglement scaling law: Any 

change of the topology will lead to a non-differentiable alter- 
ation of the prefactor of the leading order term. 

The second setting is the one of M = N/ 2 (case (i) above), 
where — in contrast to the double scaling limit — M depends 
on N, leaving one limit to consider. From the discussion 
above, we would expect £ = 00. Indeed, making use of a 
quadratic lower bound to the entanglement entropy ifl fT^l . it 
can be shown that fermions violate the area law in this setting: 
For fermionic models with nearest-neighbor interactions with 
half-filling in D = 2, so Vij = a<5dist(i ,j),i> we find after some 



straightforward algebra | 2l|] lirrijv— >oo X)&=i Es(k)/N > 
liniTv^oo Ya-Ii'^K^I) = 00. 

Bosons. - We concentrate on the geometrical setting M = 
N/2 (case (i) above) and the most significant model: The D- 
dimensional free Klein Gordon field, 



H= 1 - 
2 



dr 



[0,L]> 



7r(r) 2 + v 2 (f>(r) 2 + v 2 (V<£(r)) 5 



(4) 



which may be obtained from a Hamiltonian as in Eq. (Q~|i 
with a nearest neighbor coupling Vij = (/.t 2 + 2DQ, 2 )Sij — 
^ 2< 5dist(i,j) ,1 : Denoting the lattice spacing by a = L/N and 
taking the limit N — » 00 while keeping v 2 = fl 2 a 2 constant, 
one obtains (0J. Rescaling w.l.o.g. V 1— > V/([i 2 + 2Dfl 2 ), we 
find Vij = Sij - Civ<W(i,i),i a nd cat = (/i 2 L 2 / '(v 2 N 2 ) + 
2D) -1 jv-kxj 1/2D . Now, demanding the system to be crit- 
ical uniquely determines cm — > 1/ (2D) as the energy gap be- 
tween the ground and first excited state is given by the square 
root of the smallest eigenvalue of the coupling matrix V . 

For each individual chain k', the V(k') are transformed 
nearest-neighbor coupling matrices: 



V(k')) = \l - 2c N J2 cos(27r^/A0] Si j - c N S dmt{tj) 



In the infinite system limit, the energy gap AE((p') between 
the ground and first excited state of each chain is given by 
AE((p') 2 D = D — J2 d cos(tp' d ) — 1. Hence, in contrast to 
the critical fermionic case, the set of solutions to \ v = is of 
measure zero as only a single chain becomes critical, when we 
identify with 2ir in the spectrum, see Fig. Q] We now make 
use of a powerful result of Ref. ||4j]: The exact form of the 
logarithmic negativity E n (an upper bound to the entropy of 
entanglement and an entanglement monotone ll22i l23ll ) with 
respect to the split [1, . . . , N/2] | [N/2 + 1, . . . , N] for a har- 
monic chain with nearest neighbor coupling. It is not only an 
asymptotic statement in N, but a closed-form expression and 
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we find that limjv^oo Ylk' E n (k')/N D 1 converges to 



2(2tt) 



D-l 



[0,2tt] 



d^io g2 ( ^-g^ cos ^ +i ) 



^ - Ed=i cos (^) - 1 



independent of the mass /.t. For D = 2 it evaluates to log 2 (3 + 
2\/2)/2 and similarly for D > 2. Hence, the entanglement 
entropy for this critical model is bounded by an expression 
linear in the boundary area, we do not encounter an infrared 
divergence here, and the prefactor can be exactly determined 
for the logarithmic negativity. 

Summary. - In this work, we have clarified the issue of 
scaling of the entanglement entropy in bosonic and fermionic 
lattice systems. Our analytical argument indeed confirms and 
resolves previous numerical findings and conjectures on the 
scaling of entanglement in ground states of many-body sys- 
tems. The difference between the behavior of bosons and 
fermions may be taken as unexpected. A new length scale 
emerges that can be referred to as "entanglement thickness": 
This length scale is associated with the possibility of approx- 
imately disentangling the two regions with local operations. 
As a pure Gaussian state is equivalent to a product of two- 



mode squeezed states up to local unitary rotations, one finds 
that this scale is in fact associated with the lattice spacing, 
and not the correlation length. The violation of the area-law 
for fermions is in fact intertwined with the specific role of 
the Fermi surface. We found that quantum phase transitions 
involving an alteration of the topology of the Fermi surface 
result in a non-analytical behavior of the prefactor. 

For typical critical models we found that those individual 
boundary crossing chains obtained from our decoupling pro- 
cedure that are critical form a continuum for fermions and are 
finite in number for bosons. Integrating over all chains then 
"lifts" singularities in D > 1 for bosons. For critical bosonic 
models exhibiting only a finite number of ground states, one 
might expect that an area law holds in dimensions D > 1 even 
for models that go beyond quadratic Hamiltonians considered 
here and are thus truly interacting. Confirming or refuting this 
conjecture is an interesting challenge. 
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